Geometric phases in open tripod systems 
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We first consider stimulated Raman adibatic passages (STIRAP) in a closed four-level tripod 
system. In this case, the adiabatic eigenstates of the system acquire real geometric phases. When 
the system is open and subject to decoherence they acquire complex geometric phases that we 
determine by a Monte Carlo wave function approach. We calculate the geometric phases and the 
state evolution in the closed as well as in the open system cases and describe the deviation between 
these in terms of the phases acquired. When the system is closed, the adiabatic evolution implements 
a Hadamard gate. The open system implements an imperfect gate and hence has a fidelity below 
unity. We express this fidelity in terms of the acquired geometric phases. 

PACS numbers: 03.67.Lx,03.65.Vf,03.65.Yz 



I. INTRODUCTION 

Each cigcnstatc of a non-degenerate quantum system 
that evolves adiabatically in time acquires a well-defined 
phase in addition to the usual dynamic phase. The for- 
mer is called the geometric phase because it depends on 
the path traversed in Hilbert space [H. With the cur- 
rent interest in quantum information and computation, 
the geometric phase has received new attention because 
it is expected to be robust against some sources of deco- 
herence [1, [1,13, Hi and hence useful, e.g., for creation of 
quantum gates. To study its robustness under decoher- 
ence it is essential to generalize the concept of geometric 
phase to open quantum systems. Various proposals have 
been made 0, Q and they all point to the problem that 
phase information tends to be lost when the system is 
open and decoheres, e.g., due to spontaneous emission. 

The full system including decoherence can still be de- 
scribed, for example, by the density matrix approach, 
which predicts the relative phases between the involved 
basis states, but the information about the phases ac- 
quired by each eigenstate is not available nor is the in- 
formation about the geometric or dynamic nature of the 
phase. In this work we are interested in the phase dy- 
namics of the states of the physical system when the sys- 
tem is subject to decoherence. To this end we consider 
the Monte Carlo wave function (MCWF) approach 
We follow each wave function trajectory and calculate the 
complex geometric phases that are acquired by the adia- 
batic eigenstates. The MCWF approach has the advan- 
tage that we gain information about the evolution of the 
single trajectories and the dynamic or geometric nature 
of the phases. Not surprisingly, on average the trajecto- 
ries reproduce the density matrix result. We show how 
the evolution of the open system can be described by adi- 
abatic eigenstates that acquire complex geometric phases 
and we discuss the deviation from the closed system case 



in terms of the phases acquired. We consider a tripod 
system with three laser fields applied and present a full 
calculation of the adiabatic evolution of the wave func- 
tion of the system. The adiabatic evolution considered 
consist of a sequence of stimulated Raman adiabatic pas- 
sages (STIRAP) [loj and implements the Hadamard gate. 
We quantify the effect of decoherence in terms of the fi- 
delity of this gate and identify the role of the complex 
geometric phases. 

The paper is organized as follows. In Sec. |TT] we re- 
view the theory of geometric phases, we explain their 
appearance in STIRAP processes and present the tripod 
system. In Sec. Illll we present calculations for the closed 
system and in Sec. IIVI we derive the full solution for the 
evolution of the open system evolution and calculate the 
fidelity of the Hadamard gate. Sec. IVl concludes. 



II. BACKGROUND 

A. Geometric phase 

The adiabatic theorem states that for a given set of in- 
stantaneous eigenstates, V'n(i), and eigenenergies, En[t), 
of a timedcpendent Hamiltonian, there is no population 
transfer between the eigenstates if these vary slowly com- 
pared with the energy difference between eigenstates 
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\Er, — E„ 
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In this adiabatic regime the eigenstates do not only ac- 
quire a dynamic phase, = J^^ En{t')/h dt' , but also a 
geometric phase, 7„ 



^(0) = ^c„(0)^„(0) ^ *(t) = ^c„(0)e'(''"+^")^„(0, 
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where the geometric part of the phase can be calculated 
directly from the eigenstates [l| , 



7„ = z / {MR)\'^ iilMR)) ■ dR. 
iRi 



(3) 



Here R are time-dependent parameters of the Hamilto- 
nian, and the geometric phase becomes an integral in 
the space of these parameters (See Eq. ([7]) below for a 
specific example). 



B. STIRAP in lambda system 

STIRAP is an efficient, adiabatic process for popula- 
tion transfer in three-level systems. In the method two 
laser fields are applied to the atomic lambda system [see 
Fig. (Ha)] E [H, Q E [H, [O, p. The instan- 
taneous adiabatic dressed states for this system are two 
bright states and one dark state. The explicit form of the 
dark state {\D)) with zero energy eigenvalue, lo^ = 0, 
reads 
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no 



|2) (4) 



where and ^1-2 are Rabi frequencies (see Fig. [TJ, 
tan^ = |r2o|/|^2| and ip2 is the time-dependent phase 
difference between the two laser fields. The real ampli- 
tudes of the fields, AQ{t) and A2{t) arc time-dependent 
and the Rabi frequencies explicitly read 



n2{t) = A2{t)e-''^-^'l 



(5) 



The pulses arc modeled by sin^-pulses with different am- 
plitudes 



K,j sin ( 



2/7r(t-tsj) 



2t 



otherwise 



2r 



(6) 

where t^j is the instant of time when the pulse starts 
and T the FWHM. With all population initially in the 
|0)-state and only Q2 applied the system is in this dark 
state. Now, adiabatically increasing flo and decreasing 
^2 causes population transfer from |0) to |2), while all 
population remains in \D) of Eq. This implies that 
the excited state, |e) is never populated and hence no 
population is lost due to spontaneous emission. 

STIRAP is well-suited for studying geometric phases, 
because no dynamic phase is acquired, ui^ = 0. The 
geometric phase is calculated by Eq. ([3]) with R= {6,(p2), 



Rf 

jD = il {D\Vj^\D)-dR 
IRi 



(P2 sin^ 0dt. (7) 




FIG. 1: (a) Three-level lambda system with two laser fields 
applied with Rabi frequencies flo and ^2- (b) Pulse se- 
quence consisting of two STIRAP processes separated by AT 
in time. The first set of pulses transfers population from 
|0) to 1 2), while the second set transfers it back. We show 
the real amplitudes, ^0 and A2 of the Rabi frequencies Slo 
and 0,2 as defined in Eq. ((S)). With this pulse sequence 
sin 6 — i7o/(\/l^oP + l^aP) is zero before the first pair of 
pulses arrive, one in between the to processes and zero after 
the second pair of pulses. The FWHM of each pulse is r and 
the delay between pulses within one process is At. 



To obtain a non-zero jd, the phase difference between the 
two laser fields, (p2 , should be controlled and have a time- 
dependence with non- vanishing tp2 ■ In order to be able to 
control the actual value of we require that sin^ 9 = 
before and after the pulse sequence, ensuring that only 
during the sequence a geometric phase is acquired. As 
stated in Eq. ^ sin^ 9 = |Oo|V(l^^oP + |f^2p), which im- 
plies that during one STIRAP process sin^ 9 is increased 
from zero to one. Applying a second STIRAP process, 
where is decreased while ^,2 is increased, decreases 
sin^ 6 from one to zero. This second process transfers 
the population from |2) back to |0). The whole pulse se- 
quence is shown in Fig. [IJb) and after this sequence the 
system ends up in 



\D)=e^''-\0), 



(8) 



where 7^1 is calculated form Eq. ([7]) with = ta and 
t/ = ifc + AT as defined in Fig. \V[h). 



C. STIRAP in tripod system 

In the lambda system described above the geometric 
phase acquired by the dark state, \D) is a collective phase 
on the two atomic states, |0) and |2). We now turn to 
the tripod system shown in Fig. [H^a). This system has 
the advantage that two dark states are populated and 
these acquire different geometric phases leading to rela- 
tive phases between the atomic states as well as changes 
in their population due to changes in the interference 
of the two dark states. In the tripod system we can 
therefore calculate the geometric phase and see how it 
affects the measurable populations and relative phases 
between the atomic states. The tripod system further 
constitute a system, where a universal set of quantum 
gates can be implemented (see, e.g., [l^ and references 
therein) . The level structure consists of three lower states 
(|0), |1) and |2)) coupled to an excited state, |e), with 
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FIG. 2: (a) Four-level tripod system with three laser fields 
applied with Rabi frequencies Qq, ^i, ^2- (b) Pulse sequence 
consisting of two STIRAP processes separated by AT in time. 
The first set of pulses transfers population from {|0), |1)} to 
1 2), while the second set transfers it back. We show the real 
amplitudes, Aq, Ai and A2 of the Rabi frequencies i}o, ^1 
and Q2 as defined in Eq. The FWHM of each pulse is r 
and the delay between pulses within one process is At. 



three laser fields. These have Rabi frequencies Hq, ili 
and f22, respectively, and for simplicity we assume that 
they are all on resonance. We restrict ourselves to the 
case where the relative phase between and fli, (poi, is 
time- independent while the phase of U2, f2{t), is time- 
dependent. The real amphtudes of all fields, Ao(t), Ai{t) 
and A2{t) are time-dependent and the Rabi frequencies 
explicitly read 



rio(i) 



Ai{t)e 



(9) 



-iV2{t) 



The pulses are as in the lambda-system case modeled by 
the sin^-pulses in Eq. When we apply the pulse se- 
quence shown in Fig.[2{b), the STIRAP process transfers 
part of the population from |0) and |1) to |2) and back. 
During this process geometric phases are acquired. Cal- 
culations of the evolution of the system is presented in 



Sec. mil for a closed system and in Sec. IIVI for an open 
system. 



III. CLOSED SYSTEM 

In the rotating wave approximation, and with Rabi fre- 
quencies as defined in Eq. ([5]), we derive the Hamiltonian 



H{t) 



Ao{t) 



Ai[t)t 



A2(0e*'^^(*) 

(10) 



expressed in the {|0), |1), |e), |2)} basis. We parameterize 
the complex Rabi frequencies as 



r!o(t) - sin^oiV^oW' + ^iW', 

= cos0oiV^oW2+Ai(i)2e-'^°\ 



(11) 
(12) 

r!2(i) = cose H{t)^Aa{tY + A^{tf + A2(t)2e-'^=W, 

(13) 

where the two angles are defined as 

tan^oi = Ao(t)Mi(i), (14) 
tan0H(i) = ^Al{t)+Al{t)/A2{t). (15) 

A diagonalization of Eq. pop gives the four energy eigen- 
values 



and eigenvectors 



(i = l,2), 
(16) 



\Di) ~ — cos Oh 



sin6'H(i)(sin6'oi|0) -f cos6loie*'^oi|l)) - |e) + cos 6'//(t)e^'^2(*) |2) 



(17) 



sin 9h 



(^)(sin^^ol|o) 



-coseoie^^^Ml)) 



- cos 9h 



(t) (sin 001 1 0) + cos 001 e"^"^ 1 1) ) + sin 0ff (^e'"^^ 1 2) 



(t)e''^^(*)|2) 



17^2) = cos 001 10)- sin 0oie*'""Ml>- 



We assume that we start out in a superposition of the 
two dark eigenstates at time tj 

\DiU)) = CdAU)\Di{U)) + CD,iU)\D2{U)), (18) 

and that the evolution is adiabatic. Then the population 
stays within the space spanned by the two dark states 



and at later times the wave function is given by 

m))=CDAt)\Diit))+CDAt)\D2it)). (19) 

In order to determine the time evolution of the coeffi- 
cients {CDAt),CD2{t)} we follow Refs. Insert- 
ing Eq. p9p into the time-dependent Schrodinger equa- 
tion yields two coupled differential equations that can 
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readily be solved and we find the very simple evolution 
Cz3,(t) =e^^-iCo,(i,), (20) 

Here the phase 

ID, =- ip2sin^ Ondt', (21) 

acquired by \Di) is purely geometric because the dark 
states do not acquire any dynamic phases, = 0. 

Since no population is transferred between the two dark 
states the geometric phase could also be calculated using 
Eq. ^ . This approach was used in [l^ . 

STIRAP (Sec. Ill B[) ensures control of the geometric 
phases. The exact pulse sequence is shown in Fig. [2][b). 
The first set of pulses transfers population partially from 
|0) and |1) to |2) while the second transfers all population 
back to |0) and The amplitudes of JIq and fli arc 



such that 0Qi is kept constant. The amplitude of 

is ^max,2 = ^Aiax,0 + ^rnax,! ^ud the pulse of ^2 IS 

delayed with respect to the pulses of Qq and fli and hence 
sin 9 nit) is varied from to 1 when the first set of pulses 
arrive, while the second set of pulses adiabatically turns 
the sin6'^f(t) factor back to 0. After the whole pulse 
sequence the system ends up in the final state 

(22) 

= [- sin OoiCd, {U ) e'^^i + cos ^oi Cd, (t. )] 1 0) 
+ hcos0oiCc,(i,)e'^"i +smeoiCD2iU)]\l), 

where we have used Eq. ([T7)) in the second line. In the 
{|0), |l)}-basis, an arbitrary initial state Itpi) = ai\0) + 
bi\l) is transferred by the unitary matrix [/ to a final 

state IV'/) = Ultpi) = a/|0) + bf\l) with 



U = 



COS^ Uqi 

cos ^01 sin 9oie 



o»7-D 



1 sm 



cos Ool sm Uoie 



'Vol fg^TDi 



1) 



sin^ 6, 



01 



e 



-»V0l (g'7Di 



(23) 



By carefully adjusting the amplitudes and phases of the 
laser fields, the values of ^oii ^ai and 7^)^ can be con- 
trolled and thus generate rotations in the {|0), |l)}-basis. 
We note that U is the identity when no geometric phase 
is acquired, =0. As a special case 6*01 = ^, 0oi — 
and = —it implement a Hadamard gate 



V2 



1 1 

1 -1 



(24) 



The value 0oi = f is obtained by choosing A^, 



AmaxA^- 1), Amax,2 = \J ^l^ax .0 + ^liax ,1 ^nd ip2 = 

t/r, l\T jr = TT assures -fu-, = — tt. All simulations pre- 
sented throughout this work will use these parameters 
and the initial state = |0) which is transferred to the 
final state \ipf) = U\ip,) = (|0) + |1))/a/2 in the closed 
system case. 



IV. OPEN SYSTEM 

We use STIRAP to transfer population among |0), |1) 
and 1 2), while keeping the population in the excited state 
|e) negligible. Decoherence due to spontaneous emis- 
sion will therefore have very little effect, while dephas- 
ing caused by, e.g., collisions or phase fluctuations of the 
laser fields, will infiuence the evolution. The evolution of 
the open system can be found by solving the Lindblad 



master equation [2l| 



m m 

(25) 

where H is the Hamiltonian for the closed system and the 
decoherence is described by the Lindblad operators. Cm- 
The Lindblad master equation results in an ensemble av- 
erage of the evolution, but does not reveal a clear distinc- 
tion between the geometric and dynamic phases. We wish 
to follow the evolution of the wave functions, the acquired 
geometric and dynamic phases and how these affect the 
relative phases between and the population in the atomic 
states. Towards this end we use the quantum jump ap- 
proach, where the wave function is evolved stochastically 
H, . For a small timestcp At the evolution of the wave 
function is described as either a jump to C,n|'0(t)) or by a 
no-jump evolution with the non-Hermitian Hamiltonian 
H = H + H', where H' = -ift/2 Cl„C„. After a 
timestep with either the jump or the no-jump evolution 
the wave function is normalized. The probability for a 
jump to Cm\i>{t)) in At is P,„(t) = A<(?A(t)|Ct„C™|V(i))- 
For the method to be valid the total probability for a 
jump in At has to be small, P = J2m-^^ ^- 
method leads to many different traces, which on average 
reproduce the density matrix. For further details see, 
e.g., Refs. 

We model the dephasing by a single Lindblad operator, 
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Co = \/2ro|0)(0|, yielding the master equation 



tripod system. H = H + H' reads 



ToPoi 
To/oio 

roP2o 















(26) 



Given Co we can calculate the no jump evolution with H 
(Sec. IIV Ap as weh as the jump traces fSec lIVB]) . where 
the system is projected onto the state Cm|i/'(ij)) oc |0) at 
the instant of time t,-. 



-2iro Ao{t) 

Ai{t)e''^°^ 

Ao{t) Ai{t)e-"^o^ A2{t)e-"^^'-''> 

A2(t)e''^^'*' 

(27) 



A. Non-Hermitian no jump evolution 

In order to determine the time evolution of the system, 
With the Lindblad operator Co = V2r^|0)(0|, H' = we transform into the interaction picture with respect to 
— i/iro|0)(0| and the non-Hermitian Hamiltonian of the H' 



Ao(t)ero(*-*') 

Ai{t)e''^°' 

ylo(i)e~^«(*~*') Ai{t)e-'^°^ A2{t)e''^^'-*^ 

A2(i)e*'^2(*) 

I 



(28) 



The subscript / indicates that the evolution is described and the right (subscript r) and left (subscript /) eigen- 
in the interaction picture. This Hamiltonian is non- vectors 
Hermitian due to the Fo-exponents and in order to de- 
termine the geometric phases we follow the procedure of 
Ref. [l^]. We diagonalize Hj and find the eigenvalues 

u;f ^±yAl + Al+Aj , c.f=0 (2 = 1,2), 

(29) 



1+,)^ = -1= [sin0^^(t)(sin0oie^°^*-*'^|O) +cos0oie*'"°11)) - |e) + cos 0H(t)e^'^^(*' |2) 
V2 L 

, . 1 



1-.)/ = 



V2 



3in 6*^/ (t) (sin 6*016^° 1 0) + cos efoie**^"! 1 1)) + |e) + cos 6'H(t)e''^'^*^ |2) 
3 6'H(t)(sin6'oie^°(*-*')|0) +cos6loie*'^'"|l)) +sin6'H(i)e*'^=(*)|2), 



\Dir)l = - COS ^ 

\D2r)l = COS 0016^°^*^*'^ |0) - sin 001 



e''^°i|l). 



sin 6H{t){sm Ogie 



,-ro(t-t.)(o| +cos6loie-''^°i(l|) - (e| + cos e'H(t)e-''^^(*) (2 

= [sin0H(t)(sin0oie~^°(*~*'^(O| +cos0oie-''^°Ml|) + (e| + cos0H(i)e-''^^(*)(2| 
V2 L J 

i{Du\ = -cos0H(t)(sin0oie-^°^*-*'\O| +cos0oie""^"Ml|)+sin0H(i)e-*'^^(*)(2|, 
i{D2i\^cos9oie-^°^'-'^\0\-sm9oie-''^"'{l\. 



(30) 



The left and right eigenvectors fulfill the biorthonormal Initially (t = ti) the eigenvectors of the open system 
condition {ii\jr)i = Sij [Eq. ([50]) ] coincide with the eigenvectors of the closed 
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system [Eq. p7|) ] , and hence the initial state of the closed 
system is also the initial state of the open system 

= CDAU)\DiriU))l + CD,{U)\D2r{U))l. (31) 

The adiabatic STIRAP evolution ensures that the 
system remains within the subspace spanned by 
{\Dir{t))i,\D2r{t))i}. Inserting Eq. dST]) into the time- 
dependent Schrodinger equation gives a set of equations 
that can be solved numerically for Cdi and Cd2- With- 
out loss of generality we write the solutions as 

Cz,,(t)=e-ro"(*)e'^i(*)Cz5^(t,), (32) 
CoAt) =e-ro/'(*)e'^^(*)C^^(t,), 
and expand the wave function as 

\m)i ^e-^"''^'^e'''^Wc^^{t,)\D,rit))i (33) 

The two dark states each acquire a complex geometric 
phase composed by real (71 and 72) and imaginary parts 
{Toa{t) and ToP(t)) parameterized by the dephasing rate, 
Tq. As an example we choose an initial state = |0) 
and apply the pulse sequence in Fig. ^h) with param- 
eters leading to 9oi — |-, 001 = and -jdi ~ — ^ (see 
details in caption of Fig. [31). As discussed in Sec. Illllthis 
evolution corresponds to the Hadamard gate, Eq. pi]) 
for the closed system. In Fig. [3] we show the evolution of 
the phases as a function of time for different values of the 



dephasing rate Tq. The real geometric phases 71 (light 
grey) and 72 (grey) arc unaffected by the dephasing and 
their values are identical to the analytical result for the 
closed system [Eq. (pij) ]. which are marked with crosses 
in Fig. [3l The exponents a and /? are also almost unaf- 
fected by the dephasing rate. The results for For = 10^^ 
and For = 10~^ (solid) are identical while increasing the 
dephasing to For = 10~^ (dotted) gives only a small 
deviation. Experimentally dephasing rates can be kept 
below For = 10"'^ and for these values of Fq, a and /3 
are unaffected by the dephasing rate and hence the imag- 
inary part of the phases Foa and TqP scales linear with 
the dephasing rate. It should be noted that the value of 
a and /3 depend on the initial state while 71 and 72 are 
unaffected and equal to the values in the closed system 
case. Going back to the Schrodinger picture yields 

m)) = ^e-i^°(*-*')l")<°l[e-r«"e'^^Cz,,(t,)l^ir>/ 
V N 

(34) 

where the wave function was re-normalized (factor 
because the non-Hermitian Hamiltonian does not 
preserve the norm. In the {|0), |l)}-basis an initial state 
in the Schrodinger picture ~ ai\0) + bi\l) is trans- 
ferred to a final non-normalized state \i'f)non-norm = 
= a/|0)+5/|l), with 



L = 



COS Woie" 



-ro/3 



p*7l 



sm Woie" 



-Foa 



COS 6*01 sin6'oie*'^oi(e 



cos Oqi sm UQie 
e^''^ SIT? Oqic 



p«7i 



cos Uoie 



en2g-ro/3y 



(35) 



The normalization constant, = \CDi{U)\^e~''^^°"' + 
\CD2{ti)\'^e~'^^°'^ , depends on the initial state, and the 
final state reads \ipf) = -^lipj) non-norm- While the 
phases 71 and 72 are robust with respect to dephasing, 
the appearance of the imaginary phases, affects the pop- 
ulation of the states. Figure [3] shows the evolution of the 
population of the four states (|0), |e), |2)) with initial 
state IV'i) = |0) and the pulse sequence in Fig. [^fb) with 
parameters 6'oi = tt/S, (poi = tt and = — tt (see de- 
tailed parameters in the caption of Fig. 2]) . In the closed 
system these parameters lead to an implementation of 
the Hadamard gate and hence final populations Po = ^ 
and Pi = i. When we introduce dephasing \tpi) = |0) 
leads to final population found from the no-jump evolu- 



tion [Eq. (ESI)] 



(C0S2 6*016-^0/3 _ gjj^2 g^^g^roa-)2 

sin^ 6loie"2roa _ cos^ 6'oie-2ro/3 ' 
(sin 6*01 sin6loi(e-ro/3 + g-roa))2 
sin^ 6'oie-2roa _ ^^^2 6)g^g-2ro/3 ' 



The deviation from the closed system case (Pq = 
(cos^ 6*01 — sin^ 6*01 )'^, Pi = 4sin^ 9ni cos^ 601) is thus de- 
termined by the the two imaginary geometric phases, Foa 
and Fo/3. The values of a and /3 are found by numeri- 
cally solving the time-dependent Schrodinger equation. 
On the scale of Fig. [5] there is no deviation between the 
closed and the open system results for realistic dephasing 
rates. 
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t/x 



FIG. 3: Time evolution of 71, 72, a and /3 for different de- 
phasing rates, Fq. The calculations were made with sin'^ 
pulses (Eq. ^) and all parameters are given in units of the 

pulse width, r: (p2 = t/r, Amax.OT/'iTT = 300, Amax,l = 

A^a..o/iV2~l), = sf^l^^^Z^+AlZ^^, At/r = 1, 

and AT/t — tt. These parameters lead to 9oi = j, (poi ~ it 
and 71 = — TT. Dephasing rates For = 10~^ and For = 10^'^ 
(solid) yield identical results, while For — 10^^ (dotted) 
shows a small deviation at t/r > 3.5. 



1.0 



0.8 



.9 "-^ 

g- 0.4 



0.2 





Pi--- 
p 

e 

P„ 



-y- 



1 2 3 4 5 

t/x 

FIG. 4: Evolution of the population in the |0)-state (Po, 
solid), the |l)-state (Pi, dashed),the |e)-state (Pe, dashed- 
dotted) and the |2)-state (P2, dotted) with all population ini- 
tially in the |0)-state. The population evolve to Po = 1/2 
and Pi = 1/2 when no dephasing is present. With dephasing 
(FoTo — 10~^) the master equation as well as the no jump 
evolution [Eq. (|35}] show no deviation from the no dephas- 
ing case on the scale of this figure. The calculations were 
made with sin^ pulses [Eq. ([6])] and all parameters are given 
in units of the pulse width, r: ip2 = t/r, AmaxfiT /2-r = 300, 

Amax,l = Amax,o/{V^— 1), Amax,2 = \J -^^axfi + ^mai!, li 

At/r = 1, and AT/t — n. These parameters lead to doi = j, 
001 ~ TV and 71 = — tt. 



B. Jump evolution 

For the present choice of parameters the system jumps 
in a small part of the Monte Carlo traces to the |0)- 
state due to Co = 



If a jump occurs at tj we 
expand the wave function (|0)) in the adiabatic basis of 
the instantaneous eigenstates 



m,)) =|0) (37) 

= - COs6Hitj)sm0oi\Dir(tj)) + COa0oi\D2r{tj)) 

1 



V2' 



^9Hitj)sin0oii\+r (tj)> + I -r (<j))). 



After the jump the evolution is described by the no- 
jump non-Hermitian Hamiltonian [Eq. (|28p ]. In the adi- 
abatic basis we can describe this evolution by calculating 
the geometric and dynamic phases acquired by the four 
states. The dark part of the wave function evolves as 
described for the no-jump evolution. The bright eigen- 
states (l+r), |— r)) are separated energetically from each 
other and from the dark states, such that there is no dia- 
batic population transfer among these. The bright states 
acquire a dynamic phase as well as a complex geomet- 
ric phase that can be calculated directly Eq. ([3]). The 
complex geometric phase is the same for the two bright 
states 



(38) 



-—iVi^SHiOf^i j sin OHdt — 7: / <y52 cos Ondt 
-AVa5 + 7b 



while the dynamic phase is different for the (|+r), l~r))- 
states 



^± = t\ f ^Al + Al + Aldt'. 



The wave function at later times in the interaction pic- 
ture is 



\^{t)) =-cos6lH(ij)sin6ioie-^°"e*''i|Dir> (39) 
+ cos0oie-^°''e'^^|i?2r) 

v2 

+ J- sin 6 H{tj) sm9oie-^°^e''"'e''^~ \-r). 
v2 
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Going back to the Schrodingcr picture and calculating 
the final populations in the atomic states yields 



Po =l-\ cos^ eoie-^°^ - sin^ 6*01 cos6'H(tj)e-^°"e*'^0P, 

(40) 

Pi =1- sin^ 001 cos^ 6I01 \e-^°'^ + cos9H{tj)e-'^°°'e'^' p, 
Pg =— sin^ 6*01 sin^ eHitj)e-^^°^ sin^ 



1 

N' 
1 

iV 



P2 =^ sin^ Ooi sin^ 6l//(tj)e"2ro<5 ^,^32 



as 



where the normalization constant is given 
N = sin^ Ooicos^ eHitj)e-^^"°' + cos^ ^oie^^ro/s 
sin^ 6101 sin^ 9H{tj)e^^^°^ . Phases are only acquired from 
the latest jump time, tj. The final populations can in 
this way be calculated after each Monte Carlo trace 
and the deviations from the closed system case can be 
explained by the complex geometric phases acquired by 
the adiabatic states. When all complex geometric phases 
are carefully taken into account an average over many 
traces reproduces the numerical solution of the master 
equation [Eq. ([^5]) ]. This is shown in Fig. [51 where we 
enlarge the last part of the evolution of the populations 
shown in Fig. ID On this scale deviations between the 
closed and the open system are visible. The closed 
system (dash-dot-dot curves) yields final populations 
Po — 1/2 and Pi = 1/2 as expected for the initial state 
lipi) = |0) and the Hadamard gate applied (parameters 
are specified in the caption of Fig. O. The evolution 
of the open system is determined either by solving the 
master equation (solid curves) or by the Monte Carlo 
method. The final populations predicted by the no-jump 
evolution (dashed curves) deviates from the master 
equation (solid curves) on the order of 10^'^. An average 
over 10000 Monte Carlo traces (dash-dot-dash) reduces 
the deviation to an order of 10~^, while averaging over 
200000 traces (dotted curves) reduces it further to an 
order of 10^^. 



C. Fidelity of Hadamard gate 



The calculated Monte Carlo traces can be used to de- 
termine the fidelity of the Hadamard gate. For a given 
initial state \'>pi) we can determine the fidelity as the over- 
lap between the target (closed-system) wave function \ ipo) 
and the final Monte Carlo wave functions (no jump lipnj) 
or one jump at tj \ipj{tj))) weighed by the probability of 
each trace (no jump P„j or one jump Pj(tj)). The con- 
tributions to the fidelity from traces with more than one 
jump (Fi^mj) is negligible for realistic dephasing rates. 



.2 

a 
o 

PL, 



0.501- 


'' P 

^ V 0, j 10000 traces 


1, nj - 






' P 

( 0, p 


0.500 - 


/ Y ^ il, J 200000 traces 


P 

0, j 200000 trax:es 


0.499 - 




Pi' - 

1> P 




1, j 10000 traces 


0, nj ■ 



4.00 4.05 4.10 4.15 4.20 4.25 4.30 

t/x 

FIG. 5: Population in the |0)-state (Po, black), the |1)- 
state (Pi, grey), wlien the system is not subject to dephas- 
ing (Pi,ro=o, dot-dot-dasli) as well a when dephasing (For = 
10~^) is present. The graph shows only the final part of the 
time evolution, where the differences can be distinguished. 
With dephasing the full curve {Pi,p) shows the solution of the 
master equation [Eq. (|25p ]. the dashed curve {Pi,nj) the no- 
jump trace, while the last two curves are averages over 10000 
(fi.i 10000 tracGs dash-dot-dash) and 200000 (Pi j 200000 traces, 
dotted) jump and no-jump traces. The calculations were 
made with sin'^ pulses (Eq. (|6])) and all parameters are given 
in units of the pulse width, r: ip2 — t/r, Amax,QT/2Tr = 300, 

Amax.l = Amax,o/{V^~ 1), Amax ,2 ~ \J -^^axfi + ^maa:,l' 

At/r = 1, and AT/r = tt. 



The fidelity, accordingly reads 



(41) 



nj .non—norra. , 



^" j I (V^O IV'jinon— no7'm. (^j)) I dtj ~\'Fi rfj^j^ 

where we use the non-normalized Monte Carlo wave func- 
tions to avoid calculating the probability distributions 
directly. As a first approximation the fidelity can be cal- 
culated neglecting the jump traces 



F, 



--\{Mli^nj,non-norm)? = \ {^^Ul\L^|:,)\^ (42) 

+ 2e-ro("-+/3n.)|Cz3,(t,)nCi,,(t,)l' 

X COs(7i,„j - -l2^nj - jDi)- 



The decrease in the fidelity is governed by the geomet- 
ric phases acquired during the no-jump evolution. The 
non-normalized wave function for traces with one jump 
at the instant tj is found in three steps. The system 
evolves under the non-Hermitian Hamiltonian until tj, 
iV'jX^i)) = L\'4>j{ti)), where it acquires geometric phases 
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— Fga + Z7i and —To/? + 172- At tj the system jumps 
and \ipj{tj)) = ^/W^\0){0\'ipj{tj)). Finally the system 
evolves under the non-Hermitian Hamiltonian from tj 
to tf, \ipj{tf)) = L\'ipj{tj)), where it acquires geomet- 
ric phases —Toa' + ij'i and — rQ/3-|-i72. The fidelity from 

2 

the jump traces will hence be proportional to \/2rQ 

(43) 

where (a, 71, a', 7^, /3, 72, /?', 72) is determined by the 
geometric phases in the jump traces with jump at tj. 
Traces with more than one jump will contribute with 
terms proportional to higher orders of Tq and are there- 
fore neglected. The fidelity determined from the Monte 
Carlo traces can be compared with the Uhlmann state 
fidelity calculated from the final closed system density 
matrix po and the final open system density matrix P(ro) 

M 

Fp,^ = (Try/p^\^r^^ . (44) 

Both Fi and Fp,i give the fidelity for a given initial state. 
The average fidelity can be found by integrating over the 
surface of the Bloch sphere 

F=^J F^dn. (45) 

This averaging procedure can be simplified to only av- 
eraging over the six axial pure initial states on the 
Bloch sphere, A ={|0), |1), ^(|0) + |1)), ^(|0) - |1)), 

7l(|0)+*|l)), ^(|0)-z|l))}|2l 

In Fig. [6] we show the average fidelity as a function of the 
dephasing rate by the Monte Carlo method when only no- 
jump traces are taken into account (solid, black curve) 
and when traces with no or one jump arc included (dot- 
ted, black curve). These are compared with the full mas- 
ter equation solution (dashed, grey curve). The fidelities 
are all calculated for the Hadamard gate implemented by 
the parameters ^oi = t^/S, <y9oi = and = as in 
all previous numerical results. The fidelity decreases as 
expected when the dephasing rate increases, but is still 
acceptable when the system is subject to realistic dephas- 
ing rates. The no-jump results gives a lower bound on 
the fidelity but it is necessary to include traces with one 
jump in order to get a satisfactory accuracy for realistic 
dephasing rates. Eq. shows explicitly how the no- 
jump fidelity depends on the geometric phases acquired 
during the no-jump evolution. 

V. SUMMARY AND CONCLUSION 

We presented a method to describe the adiabatic evo- 
lution of an open system using the quantum Monte Carlo 
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FIG. 6: Fidelity as a function of the dephasing rate, Fq. The 
dashed grey curve shows the density matrix calculation, Fp 
[Eq. (|44|l ]. the dotted black curve the Monte Carlo calcula- 
tion including traces with no jump as well as one jump [F) 
and solid black includes only no-jump traces (Fnj) [Eo. I42j. 
The calculations were made with sin^ pulses [Eq. and 
all parameters are given in units of the pulse width, r: 

LP2 ~ t/T, AmaxfiT /2-K = 300, Amax,l = Amax,o/{V2 — 1), 

= \/^LZo"+^L^, At/r = 1, and AT/r = n. 



method and keeping track of all acquired phases. This 
method has the advantage that it reveals the evolution of 
the single quantum trajectories and discloses the geomet- 
ric or dynamic nature of the acquired phases in the adi- 
abatic basis. We considered a tripod system with three 
laser fields applied and calculated its instantaneous adi- 
abatic eigenstates. The tripod system is subject to a 
double STIRAP process and during this the adiabatic 
eigenstates acquire complex geometric phases. In the 
closed system case the geometric phases are purely real 
(7_Di) and all population is at all times in the two dark 
(with zero energy eigenvalues) adiabatic eigenstates. The 
acquired geometric phases create a phase difference be- 
tween the two dark states and hence generate a rotation 
in the atomic {|0), |l)}-basis. With the right parame- 
ter choice this rotation implements the Hadamard gate, 
which we have used as an example in the numerical sim- 
ulations. 

When dephasing is present we used the quantum 
Monte Carlo method, where the system either follow a 
non-Hermitian no-jump evolution during the whole time 
sequence or at one or more instants of time jump to the 
|0)-state. Before, in between and after jumps the system 
follows the non-Hermitian no-jump evolution. During 
the non-Hermitian evolution the instantaneous adiabatic 
eigenstates acquire complex geometric phases. These de- 
viate from the closed system case mainly because they 
contain non-negligible imaginary parts, which lead to a 
decay of the adiabatic eigenstates and hence influence 
the populations resulting in imperfect gate performance. 
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The specific Hadamard gate simulations show that the 
fidehty is still appreciable at realistic dephasing rates. 
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